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ABSTRACT: We use path-integral methods to derive the ground state wave func- 
tions of a number of two-dimensional fermion field theories and related systems in one— 



X 

5h I dimensional many body physics. We derive the exact wave function for the Thirring/Luttinger 



and Coset fermion models and apply our results to derive the universal behavior of the 
wave functions of the Heisenberg antiferromagnets and of the Sutherland model. We find 
explicit forms for the wave functions in the density and in the Grassmann representations. 
We show that these wave functions always have the Jastrow factorized form and calculated 
the exponent. Our results agree with the exponents derived from the Bethe Ansatz for the 
Sutherland model and the Haldane-Shastri spin chain but apply to all the systems in the 
same universality class. 



I. INTRODUCTION 



The quantization of physical systems with a finite number of degrees of freedom is usu- 
ally carried out in either the Schrodinger picture or in the Heisenberg picture. In Quantum 
Field Theory (QFT) and in Many-Body Physics the most direct approach involves the use 
of second quantization and the Heisenberg picture is natural. In this approach, the theory 
focuses on the construction of the ground state and on the computation of the Green func- 
tions. The Schrodinger picture of a QFT, and of a generic many body system, is usually 
very complicated. The wave functions, whose properties can be guessed by intuition in 
simple systems with a finite number of degrees of freedom, become very complex objects. 
In the conventional second quantized picture, the states of the Fock space are represented 
in terms of eigenstates of the number operators. In the Schrodinger approach, the states 
are wave functionals of the field configurations. These functionals are very difficult to 
define"*^ and they are very rarely used. 

In most cases, these theories are discussed in the context of some form of perturbation 
theory. The unperturbed states are usually very simple and a Feynman diagrammatic 
approach is certainly simplest in the Heisenberg picture. The situation becomes radically 
different when it comes to the study of the non-perturbative structure of these theories. To 
date, there are only two approaches to non-perturbative physics. One are the semiclassical 
and mean-field approaches. The other is the explicit construction of the exact or approx- 
imate eigenstates of the systems. While the semiclassical approach is usually presented 
in terms of path-integrals, the eigenstate approach uses the first quantized or Schrodinger 
picture. Although very rarely used in QFT, it is the conventional approach for many-body 
physics, e.g., the BCS approach to superconductivity and the Laughlin picture of the 
Fractional Quantum Hall Effect. Although much has been gained from the semiclassical 
approach^, much of what is known about non-linear theories in 1-|-1 dimensions is due to 
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the Bethe Ansatz approach"^. However, unlike the BCS and Laughhn wavefunctions, the 
Bethe Ansatz wavefunctions are so complex that their structure is almost never discussed, 
except for some very special systems^'^'^ for which the wave functions become very simple. 

In this paper we study the structure of the wave functionals of fermionic theories in 1+1 
dimensions using field-theoretic path integral methods and extract their universal content. 
We focus our attention on a number of field-theoretic model and on their application to 
the physics of one-dimensional strongly correlated systems and antiferromagnets. In par- 
ticular, we consider self-interacting theories of massless fermions, including free fermions, 
the Thirring/Luttinger model and Constrained Fermi systems (Coset models) and apply 
our results to the study of quantum Heisenberg antiferromagnets and to the Sutherland 
model. These systems have the interest since they are at a fixed point {i.e., they exhibit 
conformal invariance), their Green functions exhibit well known universal properties. 

We use the methods of references [7] and [8] to show that the ground state wave 
functionals of of the systems are determined completely by the knowledge of the Green 
functions. We introduce two representation of the wave functions: the density and the 
Grassmann representations. The density representation, which is closely related to the 
conventional particle number representation of second quantization, is given in terms of the 
density-density correlation functions. These correlation functions are exactly calculable by 
path integral methods in all the systems that we discuss which exhibit conformal and Kac- 
Moody invariance. The wave functionals have universal properties which are determined 
by the universal properties of the Green functions. The Grassmann representation gives a 
more detailed picture of the wave function in terms of fermion coherent states. 

The results presented here are directly applicable to a number of interesting models of 
strongly correlated systems in one space dimension. For instance, this result allows us to 
show that the wave functions of the Sutherland model have the same universal properties 
as those of the Thirring/Luttinger model. This is interesting since the Sutherland model 
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can be viewed as a non-relativistic regularization of the Thirring/Luttinger model and it is 
solvable by the methods of the Bethe Ansatz^. The application of the methods of the Bethe 
Ansatz to models of massless relativistic fermions, such as the Thirring/Luttinger model, 
is plagued with difficulties. Our results on Coset models are directly applicable to quantum 
antiferromagnets. It is worthwhile to stress that while the Bethe Ansatz wave functions are 
usually so complex that there structure is never studied, they too have universal features 
which can be determined directly with our methods. 

The paper is organized as follows. In section II we show how to extract the wave 
functionals from the path integral formulation of field theory. We begin with a careful 
definition of the Hilbert space. Two descriptions of the wave functionals are introduced: 
a) the second quantized or "density" picture and b) the Grassmann or Bargmann picture. 
In section III we illustrate these ideas in the context of free fermions and free bosons. These 
results had been derived previously using Schrodinger methods^'^'-*^^'^-*^. In section IV we 
apply these methods to the massless Thirring/Luttinger model. Here we give an explicit 
expression for the asymptotic behavior of the wave functional in both the density and 
Grassmann representations. Our results indicate that, while the Grassmann representation 
is a sum of terms which exhibit left- right factorization, the density representation does not. 
In addition both pictures yield wave functionals with a Jastrow form but with different 
exponents. In section V we discuss Coset models and derive the wave functionals for these 
systems. In section VI we discuss the connection between the Heisenberg and Sutherland 
models ( and their generalizations) with the Thirring/Luttinger model and show that 
the wave functions have the same long distance behavior. The connection between the 
Heisenberg model and the Coset field theories is used to derive the long distance form 
of the wave functions. We show that the asymptotic form at long distances of the wave 
function for the nearest-neighbor and the l/|rp models have the same Jastrow exponent 
(as required by universality). The conclusions are drawn in section VII. 
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II. WAVE FUNCTIONALS AND PATH INTEGRALS 



Let us begin by introducing a Hilbert space adequate to describe the states of the 
fermionic systems of interest. Our description is equally applicable to both relativistic and 
non-relativistic systems in the second quantized formulation. In either case the space of 
states is Fock space. We will introduce two representations of the wave functions: a) the 
density representation and b) the Grassmann representation. 

Let |0) be a reference state on which we will build the physical Hilbert space. In 
general we will assume that |0) is a filled fermi sea, although in non-relativistic theories is 
sometimes useful to consider the case in which |0) is the empty state {i.e., the state with 
no particles). For the case of case of non-relativistic fermions we will denote with c^{x) and 
c{x) a set of creation and annihilation operators. Then a basis for the A'"-particle subspace 
of Fock space has the form \xi, . . . ,xpf) — {xi) . . . ct(a;jv)|0). These states are eigenstates 
of the particle density operator p(x) = c^{x)c{x) with eigenvalue p{x) = X^jLi ^{x — Xj). 
We will call this the density representation and label these states by the density eigenvalue, 
i.e., |a;i, . . . , a^Tv) = \p{x))- For the case of dense systems, we will label the states by their 
density eigenvalue in the form |[p]), where p{x), the eigenvalue of p{x), is a distribution. 

Now, the density eigenstates are antisymmetric under particle exchange. However, the 
density eigenvalue is even under exchange. Thus, the density eigenvalue only specifies the 
presence or absence of a particle at a given point but says nothing about their relative 
ordering. Nevertheless, we will continue to denote the states by the density eigenvalue 
and assume that an antisymmetrization (or symmetrization, for the case of bosons) has 
been adopted. Notice that, in spite of these caveats, the states represented by the density 
eigenvalue still satisfy the exclusion principle. 

There is an alternative representation in which the antisymmetry of the states is ap- 
parent. In this representation, the coordinates which label the states are Grassmann vari- 
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ables. We will refer to this as the Grassmann representation and it is the representation 
in terms of fermion coherent states-*^-*^. Let ^i, . . . , be a set of A?" Grassmann variables. 
The coherent state is defined to be • • • > ^Af) = ^^PlZ^jLi ^jct(2^j)}|0)- These states 
are overcomplete. Let \^) be an arbitrary state of the Hilbert space. The wave func- 
tions are inner products between the states and a set of reference states. In the density 
representation the wave functions are the usual antisymmetric "orbital" wave functions 
^'(xi, . . . , xj^) = (xi, . . . , x^\'^). In the Grassmann representation, they are Grassmann 
polynomials ^'(^i, . . . ,Civ)- The orbital wave functions are the number coefficients of the 
Grassmann polynomials. 

We will now derive a set of expressions for the wave functions, for both the density and 
Grassmann representations, in terms of the correlation functions of the specific system. For 
simplicity, we will discuss only the ground state wave function. These formulas generalize 
without difficulty for excited states. For the most part they also apply almost without 
changes to both relativistic and non-relativistic systems. We will discuss first the case of 
systems of non-relativistic fermions and then generalize our results for relativistic systems. 

We begin with the density representation. Following ideas that go back to Landau's 
work on superfiuid helium^^, Dashen and Sharp (DS) have shown that density represen- 
tation follows naturally from the symmetries of the physical system. If the Hamiltonian 
has a conserved charge, such as particle number, then the density p{x) — (^{x)c{x) and 
the current j{x) = ^[c^ [x)dxc{x) — c{x)dxC^ {x)] are the time and space components of 
a locally conserved vector which obeys the continuity equation, dtp{x,t) + dxj{x,t) = 0. 
From the equal-time canonical anticommutation relations of the field operators, it is easy 
to see that the density and current operators must satisfy the equal-time commutation 
relations 

[p{x),j{x')] = -idx {S{x - x')p{x)) (2.1) 
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and 

[p{x)rp{x')] = [3i^)3{x)]=^ (2.2) 

Notice that these commutation relations are non-trivial only for dense systems, i.e., sys- 
tems for which the local density operator has non-vanishing expectation values. In other 
words, this algebra only applies to operators acting on a Hilbert space of states obtained 
by acting finitely on a reference state which represents a filled Fermi sea. We will come 
back to this point when we discuss relativistic systems. 

These commutation relations imply that the states can be either eigenstates of the 
density or eigenstates of the current but not of both. Dashen and Sharp further showed 
that p{x) and j{x) are a satisfactory set of coordinates in the sense that the set of states 
labelled by p{x) or by j{x) span the Hilbert space. In the density representation, the 
current operators are represented by 

3{x)\m) = -^p{x)^:, ( ^mp])\ (2.3) 



5p{x) 

(Recall, however, the caveats concerning antisymmetry). Since these states form a complete 
basis of the Hilbert space, it is possible to resolve the identity operator / in the form 



'^pM){[p\\ (2.4) 

up to a normalization constant (which we absorb in the measure). 

We will now follow the methods of references [7,8] to write the wave functions in the 
density representation. We begin by recalling that the absolute value squared of the ground 
state wave function in the density representation |^'o[p]p is given by 

l*o[p]P= / 'DAq exp{-i / dxi Aq{xi) p{xi)} 

r (2.5) 
^1"? w AQ\Texp{i d'^xAo{x)jo{x)}\0) 

where jo{x) = p{x). This identity can be proven by inserting the decomposition of the 
identity, eq (2.4), inside the vacuum expectation value in the integrand of eq (2.5) and 
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by carrying out the Fourier transforms. Please note that the Fourier transform is a path 
integral at a fixed time xq = 0. The operators in this expression are Heisenberg operators 
of the system in the absence of sources. The vacuum expectation value in the integrand of 
eq (2.5) can be calculated from the generating functional of density correlation functions. 
Let us denote by 2[A^] the generating functional for density and current correlation func- 
tions. 2[Aij] is obtained from the path integral of the system by coupling it minimally to 
an external, classical, gauge potential ^^,ie., 

Z[A^] = J Vi^Vi; exp{iS[i;,i(;,A^]} (2.6) 

Clearly, we have 

l*o[p]P= / 'DAo exp{-z / dxi Aoixi) p{xi)} lim Z[Ao,Ai = 0] (2.7) 

Eq (2.7) tells us that j^^oMP is determined by the generating functional of equal-time 
density correlation functions. In the next sections we will make extensive use of eq (2.7). 
This formula was derived in reference [7]. 

It is also possible to find an expression for the wave function itself, not just for the 
absolute value squared. This is done by using a generalization of the identity eq (2.5). 
In order to do that we need to be able to get information about the phase of the wave 
function, not just of its absolute value squared. The commutation relations eq (2.2) imply 
that the phase information must be related to the current operators. This leads us to 
consider the full generating function of density and current correlation functions Z[Ajj\. 
In particular, we will consider the case in which the space and time components of A^ are 
non-zero on different time surfaces and look at vacuum expectation values of the form 

Z[Ao,Ai] = (0|Texp{i J d^x Ai{x) ji{x)} exp{i j d^x Aq{x) hi^W) (2-8) 

where we must set Aq[x) — > Aq[xi)5{xq — t) first, Ai[x) —>■ Ai{xi)6{xq — t') later and, 
finally, let t' — > t. In what follows we will assume that these operations are done in the 
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order that we have just specified. This prescription amounts to a choice of order among 
the operators. By inserting now the resolution of the identity cq (2.4) inside the matrix 
element in eq (2.8) and after recalling the expression for the current operator in the density 
representation, eq (2.3), we can write eq (2.8) in the form 

j Vp exp{ J dxi dxAi{xi)p{x)j^}%[p] *o[p] exp{i J dxi Aq{xi)p{xi)} 
Vp^p + pOxAi] *o[p] exp{i / dxiAo{xi)p{xi)} 



(2.9) 

Hence, after a Fourier transform we can write 

%[p + pdxAi] *o[p] = J VAoexp{-i J dxi p{xi)Ao{xi)}Z[Ao,Ai] (2.10) 

where Z[^0)^l] is the expression calculated with the prescription of eq (2.8). We can now 
choose Ai{x) to be such that, for arbitrary p{x), p{x) + p{x)dxAi{x) = p, where p is the 
average density. We can further pick the arbitrary choice of setting ^'o[p] = 1 ( ie., a 
fixed real number). With these prescriptions we can now write an expression for the full 
wavef unction: 

*o[p] = J VAoe^-i J dxip{xi)Ao{xi)}Z[Ao,Ai] (2.11) 

provided that Ai{x) satisfies the condition p{x) + p{x)dxAi{x) = p and that Z is calculated 
in the limit (and with the order) prescribed above. 

The density representation can be applied almost without changes to relativistic fermions. 
Let us consider a relativistic Fermi field '^(a;). For the most part, we will not write the 
two Dirac components explicitly in what follows. The field ip obeys canonical equal anti- 
commutation relations. The U{1) charge current = il)'-fj_j,il), normal-ordered relative to 
the filled Dirac sea, obeys the equal-time U{1) Kac-Moody algebra 

[Jo{x), Jo{x')] = [Ji{x), Ji{x')] = [Jo{x), Ji{x')] = --dj{x - x') (2.12) 

TT 
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These commutation relations, which are the relativistic version of the non-relativistic com- 
mutation relations of eq (2.2), determine the choice of representation for the wave functions. 
Here too, we will adopt the density representation, i.e., we will pick the wave functions 
to be eigenstates of the density operator Jq{x). However, in general, the density eigen- 
value alone does not specify the state completely. For instance, for free fermions it is also 
necessary to specify the densities of left and right movers separately. In the case of inter- 
acting theories, because of the presence of the chiral anomaly, it is not generally possible 
to conserve both the right and left moving densities separately. One consequence of the 
anomaly is the existence of anomalous dimensions of various operators. In what follows we 
will assume that the regularization that we are using is such that the U{1) charge current 
is conserved even in the presence of interactions. 

We can now use the same approach that we developed above for non-relativistic sys- 
tems, and apply it to relativistic models, such as the Thirring/Luttinger model. We 
will consider the Hilbert space built on the state |0), the filled Dirac sea. The U{1) 
charge density is a diagonal operator in the basis of states x^^; . . . ; x^, a;^)} which 
denotes a set of n particles (electrons) and n antiparticles (holes), with the eigenvalue 
P(^) = Ej=l(<^(a^ - Xj) - - x'j))- The probability |*o[p]|^ that a given U{1) charge 
density profile p{x) is observed in the ground state is still given exactly by the same equa- 
tion that we derived for non-relativistic systems, eq (2.7). If we set the amplitude for the 
state with uniform (zero) density to unity, the amplitude ^o[p] is given by eq (2.11) but 
with Ai{x) given by 



This change is consequence of the fact that the Kac-Moody algebra of eq (2.12) is a little 
simpler that its non-relativistic version eq (2.2). 

We end this section with a few comments on how the Grassmann representation is 
constructed using these methods. Except for the presence of various indices ( Dirac com- 




(2.13) 
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ponents, etc. ), the method is the same for both relativistic and non-relativistic theories. 
The Grassmann representation of the wave functions can be constructed in a manner very- 
much analogous to the construction of the density representation. The main difference 
is that now we want to project the ground state onto fermion coherent state It is 

straightforward to show that this can be achieved by considering the generating function 
of fermion correlation functions Z[fj, rj], i.e., 

Z[fj,r]]= J V'4)V'4}ex.-^{iS{'4},'4}) + i J (fix[fj^|; + i/jT]]} (2.14) 

The probability for the state |[x]) to occur in the ground state is 

I*[X,X]|^ = J VfjVr]Z[fi,r]] exp{-z J dx[f]x + xv]} (2-15) 

with the same equal-time limits used in the density representation. 



III. WAVE FUNCTIONALS FOR FREE 
THEORIES IN TWO DIMENSIONS 



As a warm-up exercise, we will consider first the cases of a free relativistic scalar field 
( a theory of bosons) and free relativistic fermions, both in 1 -|- 1-dimensions. The non- 
relativistic case can also be considered without difficulty ( in any case, non-relativistic 
fermions at finite density are, in one space dimension, equivalent to relativistic fermions). 
For simplicity we will work in Euclidean space i.e., ixQ — > r. Note that besides the 
usual modifications in the partition functions iS^ Sj^j, the density-current correlation 
functions (2.1), (2.12) acquires an extra factor of i. 
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a. Free bosons 



The (Euclidean) action for a free scalar field is: 

S[<t>] = \j d^xd^(j)d^(l> (3.1) 
so that the partition function Z[J] is given by: 

Z[J] = J D(t)exp{-S[(f)] + J (fxJ.cf)) (3.2) 

or, 

Z[J] = Z[0] exp[-^ J <fx(fiyJ{x)G{x,y)J{y)] (3.3) 
with G{x, y) the 2-dimensional Laplacian Green function, 

^(a;,^) = ^log|a;-y|. (3.4) 

As explained in the precedent section, one should take for the source J: 

J = J{xi)5{xq). (3.5) 

Then, calling 

G{xi,yi) = lim G{x,y) (3.6) 
we easily find, applying formula (2.7): 

|*[(^]|2 = j DJexp[~J dxidyiJ{xi)G{xi,yi)J{yi) + i J dxiJ{xi)cp{xi)] (3.7) 



or. 



=exp[-2 J dxidyi<p{xi)G-\xi,yi)if{yi)]. (3.8) 
We have now to evaluate G~^{xi, yi) satisfying: 

j dziG-\xi,zi)G{zi,yi)^5{xi-yi). (3.9) 
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Writing: 

G~^{xi: yi) = J exp[ip(a;i - yi)] (3.10) 



one gets: 



or 



g-\p) = -2\p\ (3.11) 



G-\z) = --ij (3.12) 



SO that the ground-state wave-functional for two dimensional free-bosons is given by: 

2 4: f 1 

=exp[ dxidynp{x)- r^viy)] (3.13) 

TT J ix- yY 

which is the correct answer. 

b. Free fermions 

Let us now compare this result with that obtained for the ground-state wave functional 
of a system of two-dimensional free fermions. The (Euclidean) action for fermions coupled 
to Grassmann valued external sources fj, rj is given by: 

S[fj, ^] ~ y d^^ii^i ~ vi^ ~ i'v)- (3-14) 

(With our conventions 7;r^75 = i^.^v'-iy and 75 = ■«707l)- 
The partition function Z[f}, rf\ is then given by: 

Z%rf^ = J Di;Di:eMS[v,v]) (3-15) 

or 

Z[f], rj] = exp(- J d?x(fiyf]{x)GF{x, y)v{y)) (3.16) 
where the fermionic Green function Gp given by: 

GHx-y)--^ (3-17) 
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As before, we take rj{x) = rj{xi)S{xQ) (and the same for rj) and then we need to compute: 

GF{xi,yi)= hm GF{x-y) = -^ , . (3.18) 

yo^xo 2tt [xi — yi) 

Then, applying formula (2.15), we have: 

l*[x,x]|^ = y" Df]Dr]Z[f],r]]exp{-i J (fixixV + VX)) (3-19) 



or: 



I*[X>X]|^= / DrjDrjexpi- / dxidyir]{xi)GFixi,yi)r]{yi) 

■' (3.20) 
-i / dxi{xv + rix))- 



The quadratic path-integral in (3.20) can be easily evaluated. The answer is: 

mx,x]? = exp( J dxidyixixi)G-\xi,yi)xiyi)) (3.21) 
where Gp^{xi,yi) satisfies: 

j dziGp\xi,zi)GF{zi,yi)^S{xi-yi). (3.22) 
Fourier transforming (3.22) we get, 

Gp\z) = J ^g{p)eMiP^) (3-23) 

with 

9{p) = -271^- (3.24) 

\P\ 

We then finally have: 

|*[X,X]P = exp(-- / dxidyix{xi)^^x{yi))- (3-25) 
TT J XI -yi 

This wave functional can be written in a more familiar way if work in a subspace of the 

Hilbert space with a fixed number of particles. Indeed expanding the exponential in the 

Grassmann variables we get 

I*[X,X]P = 5^ (-^] J lUdx.dyA H {xa,{xj)Xf3j{yj)) Fo[{xj,yj}]{aA} 
n=0 " \i=l / j=l 

(3.26) 
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where j, /c = 1, . . . , n and the function Fo[xi, yi; ...]ai/3i...a„/3„ is given by 

FQ[xi,yr, ...■,Xn,yn]aiP^...ar,0r. = (7l)ai/3i-(7n)a„/J„ X det ^ (3.27) 

which is the wave function of a system of n particles and n antiparticles. 

We have to mention here that the wave function (3.25) is the same we obtain working 
canonically in the Schrodinger representation. In this approach the Dirac fermions and 
iIj are represented by differential operators in the Grassmann variables x and x acting over 
the wave functional (3.25)-'^-'^ 

'^i = ^(x^ + ho).ii^^ (3.28) 



\/2 V ^X. 

% = ^ (x. + (70).,^) . (3.29) 
They satisfy the canonical commutation relations 

{lP{x)i,i;{y)J} = iJo)^JSix-y). (3.30) 

Using the property (3.22) {i.e., G^p oc /) we can easily verify that the wave functional 
(3.25) is annihilated by the second quantized hamiltonian 



/ 



d^x^{x)^^mM = ^- (3.31) 



It might be worthwhile to reanalize the case of free fermions but considering that 
external sources couple to currents and not directly to fermion fields. What we have in mind 
is the formulation of two-dimensional fermionic models a la Sugawara, that is, in terms of 
currents = i/^JniJ so that, in our quantum mechanical language wave functionals, after 
a polarization choice, should depend say on Jq. Then, instead of action (3.14) we start 
from: 

S[Af,] = J d^xii^i ^ijj - V'TmV'^m) (3-32) 
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The corresponding partition function is, after path-integration-*^^: 

Z[A] = exp(-^ j d^x(fiyA^ix)D^^ix, y)Av{y)) (3-33) 



with: 



and 



D^r,{x -y)^ (V - V"^ du)5'^{S - y) (3.34) 



V ^{^,y) = -^\n\x -y\ (3.35) 

Following the prescription of section II we can now compute the density-representation 
ground-state wave-functional. We get: 

l*o[p]P = exp j dxidyip{xi)DQQ{xi,yi)p{yi)^ . (3.36) 

The calculation of Dqq can be done more easily in momentum space. In terms of the 
Fourier transform of the vector source Ai^{x), 

/dp ~ 
—A^{p)exp{ipx) (3.37) 

the partition function (3.33) takes the form 

Z[A] = exp(-^ J d'^pAi,{p)Di,4A,,{-p)) (3.38) 

where 

D^,M = V - ^ (3-39) 

is the Fourier transform of the current-current correlation function (3.34). We are inter- 
ested in the equal time correlation function of densities which is obtained integrating Dqq 
over the frequencies 



< P{p)p{-P) > \At=0 = [ dPO- ( 1 - 2 I 

J \ P'^+PqJ 



(3.40) 



= \P\ 
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Hence the wave function (3.36) can be written as 

|*o[p]P = exp(7r J dpp{p)^^p{-p)). (3.41) 
To return to space-time variables we use the following identity 



1 f duj 



\p\ J n uj'^ + p^ 
in order to write the Fourier anti-transform of the kernel jK as 

bl 



(3.42) 



1 



' (277)"= o;^ -I- 



.\P\. 

which is nothing but the propagator of a two dimensional massless bosonic theory. Using 
the well-known result ^ ln/i|a;| for (3.43) we can finally write for the wave function (3.41) 
in real space 

|$o[p]|2 = J dxdyp{x) In n\x-y\p{y) ^ (3.44) 

We can now find the probability of any particular density profile. For instance, we can look 
for the probability of a configuration in which the N particles are located at coordinates 
xi, . . . , a^Tv and N antiparticles (holes) are located at yi, . . . , yjq ( regardless of whether 
they are left or right movers ). For such configurations, the density profile is p{x) = 
'^^^l{S{x — Xj) — S{x — Vj))- By substituting this density profile in eq (3.44) we find the 
result 

l*0({x,};b,})P=Arnii2l^iZfzl!^lZi^^ (3.45) 

Eq (3.45) gives the probability, not the amplitude, for this configuration. We will 
show below that, up to an overall sign, the full wave function is just the square root of 
this expression. This can be verified much in the same way as we did for the Grassmann 
representation: the wave function associated with eq (3.44) is annihilated by the Hamil- 
tonian. In order to prove this statement we must write the Hamiltonian in the density 
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representation. Of course this is the Sugawara-Sommerfield form of the Hamiltonian which 
reads 

H = 27r j d'^x: p{xf + j{xf : (3.46) 

where the dots denotes normal ordering. In eq (3.46) we have restored a factor of i to the 

definition of the current j{x), which is now a manifestly hermitian operator. Looking the 
equal time commutator 

[p{x),j{y)] = -i-da:S{x-y). (3.47) 

TT 

we see that the current j{x) can be represented as a functional differential operator in the 
variable p{x) by 

j{x) = -i-d.j^y (3.48) 
Then, the Hamiltonian takes the form 

H = 2./d^. :p(.)^ + -i(a,^)': (3.49) 

The action of the current over the wave function is given by 

j{x)'ifo[p] = -i-dx I dylnp\x -y\p{y)'^o[p] 

I / 1 (3.50) 

= -i- / dy p(y)*o[p] 

■n J x-y 

showing that j{x) acts over ^q[p] as the Hilbert transform of the density. For the composite 
operator : j{x)'^ : we obtain 

: jixf : *o[p] = / dydz— ^p(y)p(z)*o[p] (3.51) 

TT^J X — y X — z 

and, using the property eq (3.22), we find 

J dx : jixf : *o[p] = -J dxpixf^oip] (3.52) 
which proves our statement 

n-^o[p]^0. (3.53) 
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IV. THE THIRRING/LUTTINGER MODEL 



In this section we will compute the ground state wave functional for an interacting 
fermionic model: the Thirring model. As we did with the free fermions we can express 
the wave functional as a functional either of the fermionic variables X: X or of the charge 
density p{x). We will start with the Grassmannian representation. 

The Lagrangian for the Thirring model coupled to external sources 77, 77 is 

C = —il^i^il^ g (i^j^ij) — fji^ — i^rj. (4.1) 

We can eliminate the quartic term in (4.1) by making a Hubbard-Stratonovich trans- 
formation. The generating functional of the model becomes 

Z[f], f]]= J DipDiPDAj^el ^{0+9A)^-\A^-fii'-i,^ <fx_ ^4 ^) 

The integral over the fermionic variables is straightforward 

Z^-n] = J DA^ det/)[^]e/^(^t^J)"'^-U^' (4.3) 

where (/^f^])"-*^ is the inverse of the Dirac operator. In order to perform the integration 
in the auxiliary field Aj^ we use the following decomposition 

Ai, = ^ {d„uj - e^,ydy^) (4.4) 

In terms of the fields (j) and lo the Green function {p[A\)~^ takes the simple form 

{P[A])-^ = e^^^+i^\~^^J^^:^e'^^'^-'^U (4.5) 

ZTT [x — yY ^ 

For the determinant of the Dirac operator we use the result (3.33)-(3.34) in terms of 
the fields and uj. The partition function (4.3) becomes 

Z% rj] = J D<t>Du;eI ^iPi'^^y'^-i^Si'PM (4.6) 
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where 

S[(P, u]= J {(1 + ^)d^4>d^,4> + ^d^^d^u} d^x. (4.7) 

(The value of the fermionic determinant is given by equations (3.33)-(3.34) up to a gauge- 
breaking term 

a I A^d^x (4.8) 



depending on the arbitrary regularization parameter a. This means that the quantization 
of the Thirring Lagrangian (4.1) leads to a one parameter family of quantum theories. Our 
election a = 0, corresponding to a gauge invariant definition of the fermionic determinant, 
is justified by the fact that this value reproduces the standard solutions of the Thirring 
model) . 

Finally if we expand the exponential term containing the Grassmann variables, we can 
perform the bosonic integrals term by term. Instead of do this is more convenient to use 
first the formula (3.3) for the square of the wave functional and leave the bosonic integral 
to the end. Hence the square of the wave functional can be written as 

I*[X,X]P = j i^0£'a;i^^i^r7e/^(^[^J)"'l^*=o)^ e-^^['^''^l ^-i S ixn+fix) (4.9) 

where 

(/)[A])~^|Ai=0 = (fib(t>{x,t)+iuj{x,t) ^ 71 ^"/50(y,t)-iuj(y,t) (4-10) 

27r (x-y) 

is the equal time propagator. 

The integration in the sources fj, r] is Gaussian and gives the result 

I*[X,X]|^ = J L>0L'a;e-/^((^[^])"'l^*=o)"'^ e-^^t-^'^l (4.11) 

(we have omitted in (4.11) a factor det[(p[A])~^|^^^o]- 1^ can be proved, using for example 
a coherent-state definition of the functional integral, that this factor is constant and does 
not play any role in our problem). 
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The kernel ((/)[^]) ^\At=o) ^ is computed identically to the one of the free fermion 
case giving the result 

((/)[A])"^Uf=o)""^ = g-i(75<?!'(a;,t)-a;(a;,t))^ 71 ^-i{^^(j)[y^t)+uj{y,t)) _ 

TT (x-y) 

Finally the integration over the fields (p and uj can be done term by term in the expansion 
of the exponential of the Grassmann variables and we have for the wave functional 

(4.13) 

where the indices a^, f3i label the Dirac components and j,k = l,...,n. The function 
-^[^1; 2/1; •••]ai/3i...a„/3„ contributions: a factor corresponding to the free fermions 

contributions (the function Fq of equations (3.26) and (3.27)) and an equal-time correlation 
function of vertex operators with the bosonic action (4.7), coming from the interaction. 

F[{xj,yj}]{aA} = ^0[{^j,2/j}]{afc/3a: 



.X 

n 



(exp( Sj(j){xj, t)) exp(^ tjCpiy^, t))) (cxp(-z ^uj{xj, t)) exp(i ^uj{yj, t))) 
j=l j=i j=i j=i 

(4.14) 

where j,k = 1, . . . ,n. The variables Sj (tj) are equal to 1 if (/3J = 1 and equal to —1 
if CKj (/3j) = 2. The correlation function (4.14) can be evaluated using the Wick theorem 
with the bosonic propagators 

(O|0(f)0(y-)|O) = - ^^^1 ,^ ln|f-y| (4.15) 

and 

{0\u;{x)u;{y)\0) = In \x - y\. (4.16) 
For example, the term F[xi,yi; x^, 2/n]i2...12 is given by 

Flxi,w;...;x„,»„]i2...i2 = IMur det ^'<' ^^'u-!t!^S " 
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where the exponent /j, is 



2 1 + g'^/TT 



(4.19) 



As we mentioned before we can also compute the wave functional of the Thirring model 
in the density representation. The first step is to obtain the generating functional of current 
correlation functions of this model. Making a Hubbard-Stratonovich transformation as we 
did in equation (4.2), the required generating functional becomes 

= j DAf,det{0 + gA-$) e"^/^' '^'^ 

where 5^ is the external source. After a shift in the integration variables — > Ajj, — ^B^ 
and using the result (3.33)-(3.34) for the fermionic determinant we get 

with Df^jy{x, y) defined in equation (3.34). We can perform the Gaussian integration in A^ 
and the partition function takes the form 

where DT^l {x, y) is the Green function of the operator 



2 

D^v{x, y) = —Df^y{x, y) + 5f^y5^{x - y) 

2 2 

= ((1 + ^) V - -dn V"^ d^)5\x - y). 

TT TT 



(4.22) 



We can compute formally 



" ^ Try . , a2X-l 



/ 2\ ~1 oo / 2\ ""^ 

^ ^ n=l 



(4.23) 
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and using the fact that that the operator d^xj ^ dy is idempotent, we get for D^^ 



D^i = ( V + ^-di^ V"' 5\x - y). (4.24) 

Finally adding the second term in the exponential of equation (4.21) we obtain 

=exp|-^-^ J d'^xd'^yBi,{x)D„r.{x,y)B^,{y)^ (4.25) 

recovering the usual result for the Thirring model: the current correlation functions are 
only modified (respect to the free theory) by the presence of a constant factor. At this 
point is straightforward to write the ground state wave functional following the same steps 
we did in the free theory. The result is 

|*o[p]P = exp|(l + ^) J dxdyp{x)liiii\x-y\p{y)^ (4.26) 

We can check this result by looking if the wave functional is annihilated by the Hamil- 
tonian. For the Thirring model the Hamiltonian in the Sugawara-Sommerfield form reads 

n = 27r(l + ^) J : p{xf + j{xf : (4.27) 

where we have restored a factor of i to the current j{x). The equal time commutator 
between the density and the current is 

[pix)Jiy)] = -ii^_a^5(a; - y). (4.28) 

and the current operator can be represented as 

Ji^) = (4.29) 

7T 9_ Op[X) 

TT 

It is easy to verify that, also in this case, the wave functional of eq (4.26) is a zero mode 
of the Hamiltonian 

7^*0 [p] = 0. (4.30) 
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As a final test of our wave functional, we can verify the vacuum expectation current al- 
gebra in the Schrodinger representation. For example for the two point correlation function 
of currents we have 

/I S 

= -^dxdy J dudv In \u - x\ In \vy\ J Dpp{u)p{v)\'^o[p\\'^ (4.31) 
1 1 1 



27r2 1 + 2! (a; - y)2 

TT 

which is the correct result. 

A final comment about the wave functionals is in order. When we tested our wave 
functionals looking if they were annihilated by the hamiltonian we assumed that they were 
real (no phases appeared when we took the square root). Our assumptions were confirmed 
when we obtained the right results, equations (3.31), (3.52) and (4.30). However , in the 
same spirit of the method developed in section I, it is possible to obtain, for a wide class 
of models, the whole ground state wave functional, i.e., with the phase factors. Let us 
explain how does this idea work for the case of the Thirring model. 

We write generically the generating functional (4.19) as 

Z[B^] = J D'il;D'il;exp V] - J d^xj^B^ (4.32) 

where '0] is the Thirring action. We can write this equation as a vacuum expectation 
value 

Z[B^] = {0\exp{- J d'^xj.Bi}exp{- J d'^xp.BQ}\0) (4.33) 
and choose the fields Bq and Bi in the form 

Bo{x,t) = Bo{x)S{t-T) and Bi{x,t) = Bi{x)6{t - t). (4.34) 

Inserting in equation (4.34) the following resolution of the identity 

I = Jdp \p){p\ (4.35) 
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we arrive to the result 



Z[B^] = J Dp(0|exp{- J (fixj.Bi}\p) exp{- J d'^xp.Bo}^o[p]. (4.36) 

Now we can use the representation (4.29) for the current in the density representation and 
we have 

exp{- J (fxj{x)Bi{x)}\p) = |p + (tt + gY'^da^Bi). (4.37) 
Hence the partition function (4.33) takes the form 

Z[B^] = J Dp%[p+ (tt + gY'^da:Bi]^o[p\exp{- J cfxp.Bo} (4.38) 

and leads to the quantity 

%[p+{7v + 9Y^d^Bi]^o[p] = J VBoZ[B^]exp{ J d^xp.Bo}. (4.39) 

This last equation allows us to read-off the phase factor of the wave functional \l/o- Using 
equations (4.21)-(4.24) we can verify that for the Thirring model the wave functional (4.26) 
is real. 



V. COSET MODELS 



As a last example of ground state wave functional of a two dimensional system, we will 
consider a Conformal Coset Model. The fermionic version of the Coset Model is obtained 
by projecting out a suitable subgroup H of the original fermionic group G. This projection 
is carried out, in the path-integral approach, introducing Lagrange multipliers which forces 
the i^-currents to vanish. The Lagrange multiplier fields can be interpreted as gauge fields. 
The total effect of the gauge fields is to constrain the fermion currents in H to zero when 
acting on the physical states^^'^^'^^. The Coset Model U{2)/SU{2) can be mapped onto 
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the Heisenberg Model. We will exploit this connection in Section VI. The generalization 
to other coset of the form U{N)/SU{n) is straightforward. 

By following the methods of the previous sections we construct the wave functions by- 
first choosing a reasonable set of labels. We will choose to work in the density represen- 
tation. Since the group U (2) has several generators, we must choose one ( or several of 
them). The constraints impose the condition that all the SU{2) generators annihilate the 
physical states. For the coset U{2)/SU{2), this leaves only one U{2) available generator, 
the "charge" 'ip^Qip. In section VI we will show that this operator is essentially the same as 
the z-component of the spin in the Heisenberg model. Thus, we are led to consider, once 
again, the generating functionals of the density correlation functions. The Lagrangian of 
the coset fermion model is given by 

2 3 
j=l a=l 

where ip-' , {j — 1, 2) are two Dirac fermions, are Lagrange multipliers and 

2 

ji = E V^Sm^V'' (5-2) 

are the SU{2) Noether currents (the a"' are the Pauli matrices). These currents annihilate 
the physical states. 

The generating functional of the correlation functions of the two diagonal densities in 
U{2) is given by 



Z[a\a2]= J Dip Dip DX exp i^- J (fxC - J d'^x{pia^ + p2a 



= J DipDipDX exp | - y d^xipj ■{id + X + a^Pi + a'^P2)ip^ 



(5.3) 



where Pi = ^4^, P2 = and = Eo=l '^m^"- We can eliminate in (5.3) the 

part of the coupling with the sources by a shift in the variable A^. Then the partition 
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function (5.3) can be written as 

Z[a\ a^] = J DXdet (^^■{id + X + ^—^i^ 

f 1 /■ ,2 ,2 a'^(x) + a^(x) ^ , ,aVy) + a2(y)1 ^ ^ ^ ^^^^^^^ 
= exp|-- y ^^00(a;,y)^^^^ ^1 y exp{-Vr[A]} 

(5.4) 

where -Doo(2^) J/) is given in equation (3.34) and 1^[A] is the Wess-Zumino-Witten action 
Eq (5.4) shows that the non-abehan structure factorizes and decouples. This is a conse- 
quence of the simple structure of this coset. In more complicated coset, such as those in 
the sequence of minimal models, the non-abelian structure does not decouple. 

The factorized structure of eq (5.4) has important consequences for the wave functions. 
Indeed, once the generating function is known, the methods of the past sections tell us 
that the wave functional of the U (2) / SV (2) Coset Model is 

|*[P1,P2]P = y Da'Da^ Zla" ,a^\ ^^P dx {pia^ + p2a'^)^ 

/■ ^ 1 ^ 2 f 1 f , , + ^ a^ + a'^] 
= / Da Da exp < / axay Dqq > 

( f \ + a^ — a'^ 

exp < / dx (pi + p2) — ^ \- (Pl - P2) — ^ — 



(5.5) 



= S[pi - p2] exp |2 j dxdy pi{x) In \x - y\pi{y) 

The factor S[pi — P2] is a result of the constraint and it simply means that the two densities 
have to be equal to each other. The second factor implies that the wave function of 
the coset, as a functional of one of the densities (say pi) is equal to the square of the 
wave function for free fermions ( see eq (3.44)). Thus, the Jastrow exponent u for the 
U{2)/SU{2) Coset Model is v = 2. In section VI we will make use of this result in the 
context of the spin chains. 

It is straightforward to generalize the result of eq (5.5) for the slightly more general 
case of the coset U{N)/SU{n). In this case, the constraint affects only n — 1 diagonal 
generators of U{N) {i.e., those which span the Cartan subalgebra of SU{n)). If we define 
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an array of densities pi, . . . , p]\i, the constraint forces the first n densities to be equal to 
each other and, say equal to pi, while the other N — n densities remain unconstrained. 
The wave function ^^qIpI) • • • > PTv] becomes 

|*0[P1,---,P7V]|^ = |*o[pl,---,Pl,Pn+l,---P7v]P 

= exp |n j dxdy pi{x)hi\x — y\p\{y) 



(5.6) 



JJ ex.^ I j dxdy pi{x)\D.\x - y\pi{y)\ 



l- 

From this results we conclude that for all the simple coset the Jastrow exponents are 
always integers. Non-integer exponents can only arise as the consequence of the non~ 
abelian structure. 



VI. APPLICATION TO ONE-DIMENSIONAL STRONGLY 
CORRELATED FERMI SYSTEMS AND ANTIFERROMAGNETS 



We will now discuss the applicability of the results of previous sections to a number 
of models of one-dimensional many-body physics. A number of these models are exactly 
solvable through the Bethe Ansatz. These systems are said to be integrable and for this to 
happen a number of conditions have to be satisfied. While these conditions do not affect 
the low energy behavior, they do change the microscopic properties of these systems. The 
Bethe Ansatz wave functions thus appear to depend on a number of microscopic properties 
which are physically irrelevant but are essential for integrability to hold. We will show in 
this section that the wave functions of these systems also have universal features. These 
universal properties, independent as they are of the short distance physics, hold separately 
from the integrability conditions. Most of the identifications of models listed below have 
been known for some time. We present them here with the sole purpose of making the 
connections clear. 
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There is a considerable number of systems in one-dimensional many-body physics 
whose long distance behavior is governed by the Thirring/Luttinger model, i.e., they 
are in the same universality class. The reason behind this universality has been known for 
a long time 19'20,21. gp^^^ exchange is ignored, the low energy physics of all these systems 
reduces to the dynamics of fermion excitations with momentum near the two Fermi points 
(right and left movers) interacting via backscattering processes ( i.e., processes which 
exchange left and right movers). The basic assumption here is that, if the interactions 
are not too strong, the only states which participate in the dynamics are those which are 
sufficiently close to the Fermi "surface" . With this assumptions, it is natural to separate 
the electron field into its fast and slow components. At weak coupling, the fast components 
can be neglected (their only effect is to renormalize the coupling constant and the energy 
scale for the slow components ). If we denote by ip{x) the "electron" field operator, the 
left iI^l{x) and right i1^r{x) moving components of the Fermi field are defined through the 
decomposition 

^(x) = e^^^^ ^r{x) + e-^^^^ (6.1) 

and can be arranged in the form of a two-component spinor i/jai^) (o^ = 1) 2), with i/ji = i/jji 
and V'2 = V'L- At low energies the non- inter acting dispersion curve can be linearized in 
the vicinity of the Fermi points. The result is a spinor Fermi field with an effective low 
energy Hamiltonian which is identical to the Hamiltonian of the relativistically invariant 
Thirring/Luttinger model. This identification holds provided that (a) all energies are 
rescaled so that the Fermi velocity is set to be one , vp = 1, and (b) all irrelevant operators 
are neglected {i.e., operators with scaling dimensions larger than 2). 

a) The Anisotropic Antiferromagnetic (XXZ) Heisenberg spin chain: 

The best known example is the anisotropic Heisenberg antiferromagnet (the XXZ 
model) as a function of anisotropy. The physics of the Thirring/Luttinger model coincides 
with the long distance physics of the XXZ model-*^^ if umklapp processes can be neglected 
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In this model spin-^ degrees of freedom are arranged on a linear chain with N sites. 
The Hamiltonian is 

N 

H = jJ2 (^l(i) ^lO' + 1) + ^2(j) (r2ij + 1) + A asU) a^U + 1)) (6.2) 
i=i 

where {crj^} {k = 1, 2, 3) are the three Pauli matrices, J is the exchange constant and A mea- 
sures the anisotropy (A = 1 is the Heisenbcrg antiferromagnct). This model is integrable 
through the Bethe Ansatz 23,24 rpj^^ identification with the Thirring/Luttinger model 
was done by Luther and Peschel-*^^ (a detailed argument is reviewed in references[20,25]). 
The mapping is done in two steps. First, a Jordan- Wigner transformation is used to map 
the spin model to a model of spinless fermions c{j) defined on the same one-dimensional 
lattice. The Hamiltonian for the equivalent fermion model turns out to be 




where n{j) = c'{j)c{j) is the fermion occupation number, and it is related with the z- 
component of the spin through <J-^{j) = 2n{j) — 1. Hence, the state with zero (one) 
fermion at j is the same as the state with one spin down (up) at j. The sector with Sz = 
is identified with the half-filled sector of the fermion model. The continuum limit of this 
fermion model in this sector is easily seen to be equivalent to the Thirring/Luttinger model 
once the Fermi field is split into its left and right components. For small A {i.e., close to 
the Xy-model limit), the coupling constant for the g Thirring model is seen to be given 
hy = X (in units such that vp = 1). The relation between the two coupling constants 
is not universal and, in general, it is non-linear The umkalpp operators are irrelevant 
close to the XY limit ( small A) but become marginal at A = 1, the isotropic point. In the 
Thirring/Luttinger representation these operators become marginal at g^ — tt. For g^ > it 
these operators become relevant and cause the system to become massive. 

Since the charge density of the Thirring model is proportional to the average of the 
fermion occupation numbers on two consecutive sites {i.e., the average of the z-components 
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of the spins), the wave function of the Thirring/Luttinger model in the density representa- 
tion has to coincide with the wave function of the Heisenberg chain in the representation. 
This is the conventional representation of the Bethe Ansatz wave functions in which the 
coordinates represent spin flips. On the other hand, the Grassmann representation of the 
Thirring wave functions have coordinates which represent fermions, not spin flips. The 
elementary excitations of the Heisenberg chain are fermions which, in the spin picture, 
are solitons. In section IV we found that the wave functions obtained through these two 
pictures yield different exponents. With this identiflcations, we can now make use of our 
results from section IV and conclude that, for large separation of its coordinates, the wave 
function of the XX Z model also has a Jastrow form with the exponent that we found for 
the Thirring model. This exponent is accurate only close to the XY limit. The behavior 
of the wave function at the isotropic Heisenberg limit can be determined as a limit — > tt. 
The exponent of the Jastrow factor of the wave function, u = {1 + g'^/ir), has the limiting 
value v = 2. Interestingly enough this is precisely the same exponent found by Haldane 
and Shastri in the spin chain. 

b) Sutherland's model: 

Sutherland's model is a system of non-relativistic spinless fermions, at flxed density, 
interacting via a pair interaction potential V{\x — x'\) = A/\x — x'\'^, where A is a coupling 
constant. In second quantized form the Hamiltonian for the Fermi fleld i(^{x) is 

H = jdx ±^{d^^^(x)) (d^iPix)) + ljdxj dx'pix)Vi\x - x'\)pix') (6.4) 

where p{x) = tjj^ {x)tjj{x) is the density operator and we have set ^ = 1 and M is the mass. 

By means of a suitable generalization of Bethe's Ansatz, Sutherland^ was able to flnd 
the complete spectrum of this model. The wave functions turned out to have a very simple 
form. For a system with periodic boundary conditions it is natural to think of the particles 
as being on a ring of circumference L. For the pair potential V^(|a;|) = 
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sm(-l--L) 



which has the required periodicity properties, he found that the wave function for the 
ground state in the AT-particle sector has the Jastrow form 

^ = ll[ijix,-Xj)]^ (6.5) 
i>j 



with 



1 + (1 + 2^)2 (6.6) 



We will use the methods of the previous sections to calculate the asymptotic long 
distance behavior of the wave function. Up to a redefinition of the coupling constant, we 
will be able to reproduce Sutherland's result, eq (6.6). 

For a general pair potential V{\x — x'\), Sutherland's model is equivalent (at low ener- 
gies) to a Thirring-like model with Hamiltonian 

Iff,, , > , '•''■'''> 

- dx dx {ip' {x)'ip{x))u{\x — X \) {ip^ {x )'ip{x )) 



2 „ 

where have separated in the Hamiltonian the interactions into forward scattering backward 
scattering processes. The forward scattering processes are weighed by the low momentum 
components of the pair potential (low compared with 2kp) ^(l^^l) = ^d^^D- The back- 
ward scattering processes, which involve momentum transfers near 2kp, are represented 
by the Thirring term with the coupling constant ^ —■^V{2kp), where vp is the Fermi 
velocity and V{2kp) is the Fourier transform of the pair potential at 2kp. This separation 
is accurate for processes in which the single particle energies differences are much smaller 
than an energy cutoff D in which range the dispersion curve is accurately linear. This 
restricts the validity of this identification to the weak coupling regime. For larger values of 
the coupling constant, it is necessary to renormalize the parameters of the Thirring model. 

With these caveats we can, in principle, use the wave functions of section IV, and use 
these results in the long distance limit of the Sutherland model. But, before we do that, 
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we must find out if the density-density coupfing could possible change the results. There 
are several ways to see that the pair potential cannot affect the long distance behavior. 
Firstly, by direct bosonization, it is easy to show that, for potentials u{\x\) which whose 
Fourier transform u{q) satisfies q u{q) ^ ( as — > 0), the effective bosonized action 
coincides with the bosonized action of the Thirring model, up to a finite renormalization 
of the Thirring coupling constant (i.e., the backscattering amplitude) . Instead of following 
that approach, we can compute the wave function for the generalized theory of eq (6.7) 
using the same methods as in the previous sections. An explicit calculation in the density 
representation yields the wave function for the ground state to be 




where p{p) is the Fourier transform of the density. This is the exact wave function for this 
model (within the approximation of a linearized dispersion curve) . By direct inspection of 
this result it is apparent that the Jastrow exponent A is determined by the low momentum 
behavior of the potential u{q). Any delta function contribution to wda;!) is a constant 
term in the Fourier transform at all momenta. Thus, delta function pieces shift {i.e., 
renormalize) both the Thirring coupling or backscattering amplitude (oc V{2kp)) and the 
— > behavior of u{q), the forward scattering amplitude. Thus, such effects are essentially 
unobservable since these effects cancel each other out. If the behavior u{q) — > const, is 
excluded, we have to consider a polynomial {i.e., analytic) behavior, such as u{q) ~ A |g|. 
This is in fact the behavior for Sutherland's model where the constant is ^ = —ttq (for 
F(|a;|) = ga/\x\'^ and a a short distance cutoff). 

Obviously, this analysis fails if the Fourier transform of the pair potential is singular 
in the q ^ ^ limit in the form u{q) A jcf (with A < 0). For this type of behavior, the 
kernel in the integral of eq (6.8) behaves, at large separations, like — const. . Thus, 

while the wave function retains its Jastrow form, the physics is completely different: for 
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all potentials with cr > {i.e., they decay more slowly than at long distances), the 

wave function vanishes like exp(— const. |a;|'^/^) at long distances. This indicates that, for 
these potentials, the spectrum of these systems does not contain particle-like excitations. 
Indeed, the case a = 2 (i.e., u = —e^/q^ corresponds to the one-dimensional Coulomb 
interaction This case is exactly equivalent to electrodynamics in one space dimension, the 
Schwinger model. The Schwinger model is known to be a confining theory and it does not 
contain any states in its spectrum with fermionic quantum numbers. 

Thus, we conclude that, for a large class of pair potentials, the wave function at long 
distances not only has the Jastrow form, but the Jastrow exponent can be calculated 
directly from the Thirring/Luttinger model. In section IV we got a wave function with 
the same Jastrow form but with the exponent Aj^/j = 1 -|- g'^/ir (here g is the coupling 
constant of the Thirring model, i.e., g ^ —■2^V{2kp)). In the Sutherland model, the 
Fourier transform of the pair potential V^(|a;|) = g/\x\^ is t^(|p|) = —{g/vp)\p\. Hence, we 
must identify the Thirring model g'^/ir with Sutherland's coupling constant g (for vp = 1). 
This identification is valid for g small and, in this limit, the exponent of the wave function 
calculated in the Thirring model agrees with Sutherland's exponent eq (6.6). However, 
what matters here is that the form of the wave function is the same. The actual value of 
the exponent depends on the identification of the coupling constants, which is not universal. 
This problem is well known in the context of the Heisenberg spin chain. 

c) The Haldane-Shastri model: 

The Haldane-Shastri (HS) model is a non-local modification of the Heisenberg spin 

chain with Hamiltonian 

N 

HhS = E "^(1^' - ^1) ^l(^) + ^2(j>2(fc) + A a3(j) as{k)) (6.9) 

j,k=l 

with J(|j — A;|) = A/\j — fcp (for j ^ k). It has the same physics as the Heisenberg spin 
chain except that the Bethe Ansatz has a much simpler structure: the wave functions have 
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a Jastrow form. But, since the interactions are non-focal, the mapping of the HS model 
to the Thirring/Luttinger model cannot be done directly, contrary to what we discussed 
for the Heisenberg spin chain. 

The HS model can be viewed as a strong coupling limit of a generalization of a Hubbard 

model. Consider a system of spin-^ fermions on a one-dimensional lattice of length N with 

creation and annihilation operators cj-(j) and Ca(j) defined at each site j (= 1, . . . , N) and 

for each spin orientation a =t, j. The Hamiltonian 

TV N 
H=J2Y1 [t{\j-k\)cUj)ca{k)+h.c]+J2Un^U)HU) (6-10) 

is a Hubbard model with a non-local hopping term with amplitude t{\j — k\). Standard 
arguments^^ show that, for a half-filled system and in the strong coupling limit, (U ^ 
maxj — A;|)), this Hubbard model is also equivalent to a Heisenberg antiferromagnet. 
The exchange interactions are now also non-local and have an exchange constant 
Jilj ~ k\) = [t(\j — k\)]'^/U. Hence, the J/\x\'^ interaction of the HS model requires that 
the equivalent Hubbard model should have a hopping term t{\x\) a; yJU /\x\. Hence, our 
results on wave functions also apply to the long distance behavior of the spin component 
of the wave functions of the Hubbard models. 

However, all quantum antiferromagnets in all dimensions and for all lattices, can also 
be regarded as the strong coupling limit of a gauge theory^^. For a one-dimensional 
chain with exchange constant J{\xj — xj^]) the correspondence goes as follows: a set of 
fermion degrees of freedom are defined on each site of the chain with the corresponding 
fermion creation and annihilation operators, ipa{j) and ipa{i), a — 1,2. The spin operators 
are defined to be S+{j) = ^eabMmijl S-{j) = S+{j)^ and S^ij) = M)Ha{j) - 1- 
These operators are invariant under the local SU (2) gauge transformation V'a(j) ~^ i^aU) — 
UabU)'4'b{j) ■ terms of these fermions, the spin operators have a local symmetry, the 
quantum antiferromagnets must be equivalent to an SU{2) gauge theory. Indeed, it is 
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straightforward to show that they are the strong couphng hmit of an SU{2) lattice gauge 
theory^^'^^. The gauge degrees of freedom are 2x2 SU{2) matrices. In the absence of 
gauge fields, the free fermions have a global U (2) invariance. The effect of gauge fields is 
to set all SU (2) charges and currents to zero on the physical state^^. Hence, in the (naive) 
continuum limit, this is a theory of constrained (two-component) Dirac fermions in 1 + 1 
dimensions with the effective Lagrangian C — i/jipilj, where Djj, is the covariant derivative 
in the fundamental (spinor) representation of local SU{2). This result can be derived by 
a straightforward application of the methods used above. The constraint implies that the 
physical degrees of freedom live in the coset G = U{2)/SU{2). (We have assumed that 
any non-locality is sufficiently weak so that the mapping holds; it certainly works for the 
usual nearest neighbor Heisenberg model and it applies to the HS model as well). More 
complex coset models are found for higher spin representations as well as for groups other 
than SU{2) (for a review on quantum chains see reference [28]). 

Since all the local SU{2) currents and densities annihilate the physical states, it is 

obvious that they cannot be used to label the states. Thus, the labels of the wave functions 

have to be the eigenvalues of the unconstrained currents. The coset G — U{2)/SU{2) 

contains a global U{1) subgroup whose generator is, precisely, i/^^V which, as we discussed 

above, is essentially the same as the spin operator S'3 in the Heisenberg model. Hence, 

the wave function of the constrained fermion, on the diagonal generators of the coset 

G — U{2)/SU{2), has precisely the same information as the Bethe wave functions for the 

Heisenberg model. The results of section V show that these wave functions also have the 

Jastrow-like factorized form and that the exponent is also equal to 2 ( see eq (5.5)), in 

agreement with the results of Haldane and Shastri. 
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VII. CONCLUSIONS 



In this paper we have used path-integral methods to derive the ground state wave 
functions of a number of two-dimensional fermion field theories and related systems in one- 
dimensional many body physics. We derived the exact wave function for the Thirring/Luttinger 
and Coset fermion models and used these results to derive the universal behavior of the 
wave functions of the Heisenbergmodels and of the Sutherland model. We found explicit 
forms for the wave functions in the density and in the Grassmann representations. For the 
field-theoretic models, the wave functions can be calculated exactly. We showed that these 
wave functions always have the Jastrow factorized form and calculated the exponent. The 
exponents obtained with our methods agree with the exponents derived from the Bethe 
Ansatz for the Sutherland model and the Haldane-Shastri spin chain. We also showed 
that they also apply to all the systems which are in the same universality class. Thus, 
the long distance form of the ground state wave function for the nearest neighbor Heisen- 
berg antiferromagnet has the same exponent u = 2 as the HS model. We also verified 
that this result also applies to the fully SU{2) invariant U{2)/SU{2) constrained (Coset) 
fermion field theory. Also, up to a finite renormalization of the coupling constant, the 
Thirring/Luttinger form of the wave function agrees with the wave functions of both the 
Sutherland model and of the anisotropic (XXZ) antiferromagnet. 

Our results confirm that the wave functions of these theories also have a universal 
content. In previous studies, the universality was always investigated in the context of 
the Green functions. Naturally, the universality of the Green function requires that the 
wave function should also have universal component. However, due to the complexity of 
the wave functions, their universality is not usually investigated. Our methods can be 
extended with only minor effort to the computation of the wave functions of the excited 
states. More interesting, however, is the study of the properties of the wave functions to 
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more complex theories with Kac-Moody level A; > 1. This work is in progress. 
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